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Abstract 

^ • We obtain the contributions to the renormahzation group functions 

O \ of all the diagrams containing the unique one-loop primitive divergence 

of a simple supersymmetric Wess-Zumino model, up to more than 200 
loops. The asymptotic behavior of the coefhcients in the expansion of the 
CNJ ■ anomalous dimension is analysed. 

> 

^ ■ 1 Introduction 

O 

Perturbative Quantum Field Theory is known for its tremendous successes, with 
■ its ability to obtain highly precise values in Quantum Electrodynamics, or to 

OO , test the Standard Model with radiative corrections to weak interactions. How- 

ever, actual calculations become rapidly cumbersome and display a conjunction 
of analytical and combinatorial difficulties. The fast growth of the number of 
^ . relevant terms gets compounded by the need to subtract a growing number of 



o 



H ' sub divergences. 



In the quest of organizing principles for taming the combinatorial problem, 
a major progress has been the recognition of a Hopf algebra structure in the 
renormalization of Quantum Field Theories [U El EJ Hj. The cohomology of 
the introduced Hopf algebras has been related to Schwinger-Dyson equations 
(see [7] and the references therein). An early application to the summation 
of a category of diagrams in a simple Yukawa field theory or a 0^ theory in 
6 dimensions has been obtained by Broadhurst and Kreimer in [5]. It was 
however later shown that in this simple case, the Schwinger-Dyson equation, 
which is linear, could be solved exactly, without any reference to the Hopf 
algebra of diagrams, since it gives rise to a non-linear differential equation for 
the anomalous dimension [6]. 
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In a preceding work we adapted this work to a supersymmetric model 
with a simple Wess-Zumino type interaction. As it is well-known, in view of its 
ability to limit renormalization phenomena, supersymmetry has allowed to re- 
solve a number of theoretical problems. Although its relevance as a property of 
the physical world is up to now as elusive as ever, it is still a leading candidate 
for physics beyond the standard model. Of course, if supersymmetry plays a role 
in particle physics, it has to be broken through some (preferably spontaneous 
or dynamical) mechanism. In this respect, theories with metastable supersym- 
metry breaking have been recently very actively investigated, particularly some 
generalized Wess-Zumino models [9l [10] . 

The method employed in [6], if really powerful, does not generalize to more 
complex situations. A strategy has been proposed in ^] to deal with non- 
linear Schwinger-Dyson, in which the formulation of the renormalization group 
in terms of the Hopf algebra of diagrams becomes essential. We here apply 
this proposal in a supersymmetric setting. With respect to the cases envisaged 
in [11], the supersymmetric case is simpler because we need only consider a 
unique anomalous dimension, which is common to all the member of the su- 
permultiplet. We further simplify the procedure by a direct computation of 
the renormalized propagator, so that we can obtain explicit results up to high 
orders. 

We will first review the bases of this procedure, which are the following 
two results: (i) the renormalization group is a one-parameter subgroup of the 
group of evaluation functions on the Hopf algebra of graphs [3] and (ii) some 
combinations of graphs form sub-Hopf algebras. We will show that these results 
allow to recover the full propagators from the renormalization group functions. 

We then apply this calculation scheme to the Wess-Zumino supersymmetric 
model, building on our preceding work [8]. As we shall see, the efficiency of the 
method allows to reach high degrees of perturbation theory, being essentially 
limited by the size of the expression in terms of values of the Riemann func- 
tion. Switching to numerical approximations allows to go further and reach the 
asymptotic behavior of the coefficients of the development of the anomalous 
dimension 7 in powers of the "fine coupling constant" o. This behavior allows 
for the definition of a convergent Borel transform of this serie. A singularity is 
however present on the positive axis but we discuss the options for nevertheless 
obtaining a sensible resummation. This work suggests further developments 
that we explore as a conclusion. 

2 Exponentiation of the renormalization group 

Our purpose is to reach high order of perturbation theory in the evaluation 
of renormalization group function. Along the lines suggested in [11], we use 
Schwinger-Dyson equations to produce a family of diagrams for which analytical 
evaluation is very simple. However, a naive approach to this evaluation involves 
a triple serie in the coupling constant g, the logarithm of the momentum L = 
log((7^/^^) and the dimensional regularization parameter e. It was shown in [6] 
how to reduce the problem to a single non-linear differential equation for the 
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renormalization group 7-function, but the proposal relied on the linearity of the 
Schwinger-Dyson equation. 

Now, if we work only with renormalized quantities, we can avoid the e- 
parameter. The momentum dependence can then be recovered from the renor- 
malization group. Consider the ratio of the full propagator with the free propa- 
gator which, in a Lorentz invariant scheme and for a massless theory, is a single 
number. It can be fixed to 1 at the scale /i^ in a "physical" renormalization 
scheme and the value for any momentum p can be inferred from the action of 
the renormalization group with parameter t = \og{p'^ / y?) . Other renormaliza- 
tion schemes like the minimal substraction one or its variants would introduce 
finite contributions at the reference momentum. Similarly, vertex functions, 
which generically depend on multiple scales, would involve finite contributions 
at the reference scale. However, here we are interested only in the propagator 
of a massless theory where these complications are not present and the full 
propagator can be recovered from the renormalization group. 

In it was shown that the renormalization group is a one parameter 
subgroup {Ft\ of the group of characters of the Hopf algebra of graphs, which 
means: 

Ft{XY) = Ft{X)Ft{Y) 

Fs+t{X) = {Fs(^Ft){AX) (1) 

Multiple derivations of this identity allows to recursively obtain the derivatives 
of Ft at the origin. Indeed: 

9!'^^Ft\^^Q = d^dt^F^^^l^^Q^^Q = ^sF^I^^q S^'^FJ^^q (2) 

Applied on a group-like element of the Hopf algebra, that is, on an element X 
such that AX = X X, this identity shows that the successive derivatives are 
powers of the first one, so that Ft{X) = e"*. We will deal however with more 
general elements of the Hopf algebra of diagrams, so that the results will be 
more complicated. 

Finally since Ft is a character, its evaluation on a product is the product 
of the evaluations, so that /? = dtFt\t=o, the generator of the renormalization 
group, is a derivation. This will be important for the applications of eq. ([2]). 

3 Hopf algebra of Green functions 

In order to make the best use of the Hopf algebra technique, it is convenient to 
work not at the level of the individual diagrams, but on some sums of diagrams. 
If the sums of interest can be shown to form sub Hopf algebras, this allows to 
overcome in part the combinatorial complexities of calculations. 

Our starting point is that the effective coupling is a map from the Hopf 
algebra of the graphs to the Hopf algebra of formal diffeomorphism tangent to 
the identity in 0. For couplings associated to a three-point vertex, it will be 
convenient to express this fact not in terms of the coupling g, but of a "fine 
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structure constant" a = Cg^, in order to avoid square roots. The proportion- 
ality constant C can be choosen as in QED to simplify further results. If the 
vertex couples particle species j, k and I, the effective constant is defined as: 

Oeff = (3) 
i 2,j-L 2^ 2,1 

Here T2J is the two-point function for the j specie while r„ is the effective 
interaction vertex. 
Writing 

Oeff = ^ Oeff.na" (4) 
n 

the action of the coproduct on Oefr is given by [U [l2] 

AOeflF = OefT Oeff.n (5) 
n 

This corresponds to the fact that a^s is a map to the Hopf algebra of (coordi- 
nates for) formal diffeomorphisms. 

In fact, two refinements are easy to obtain. First, the same formula holds 
whenever we do not take for the F's the full perturbative expansion, but only 
the terms which contain a certain set of primitive divergences. This holds 
in particular if we take for the F's the approximation resulting from a given 
truncation of the Schwinger-Dyson equations. 

An examination of the proof of the previous identities in [12j shows that 
stronger results hold for the two- or three-point functions which are involved, 

^2,m = ^r^.Tifl", (6) 

n 

r. = 5^r,,„a", (7) 

n=0 

AF2,m = ^ ^2,macS ® ^m,n (8) 
n 

AF„ = ^F„a^ff®F^,„ (9) 

n 

In the cases in which the coupling constant is not renormalized, Oeff = a, these 
equations indicate that the F's are group like elements of the Hopf algebra. 
Now, when deriving relation ([5]) from relations ()8|9p . one uses the fact that 
formulas of this kind are multiplicative. Indeed, suppose that A = AnoJ^ and 
B = BnoP" are two series satisfying a property like ([8]): 

A^ = ^^<ff(g)^„. (10) 

n 

Then the product AB satisfies the same property. We can further remark that 
the coproducts of the propagators satisfy similar equations. This easily follows 
from the fact that the identity element of the Hopf algebra trivially satisfies (jlOp . 
This allows to prove ([5]) from eqs. 
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4 Practical formulas 



Joining the results of the two preceding sections allows to construct the full 
renormalized propagator at any order in perturbation theory from the first 
derivative with respect to L, that is, the renormalization group function 7. In 
fact, we can get at will the propagator or its inverse, the IPI two-point function. 

Let us consider the evaluation of the renormalization group on the propa- 
gator r^j, normalized by the free propagator, 

00 

Ft{r-^.) = i + Y,ikt' (11) 

n=l 

where, as before, t = log(p^/^^). Notice that the coefficient 71 in expansion 
pip is related to the anomalous dimension 7, 

71 = -27 (12) 

The successive terms of the development of this expression in powers of t can 
be computed from eq. ([2]), using the coproduct in eq. ([8]), 

{k + l)7fc+i = Y^hi + n^)a"7fc,„ = (71 +(3adahk (13) 

n 

To obtain the preceding equation we made use of the simple differentiation rule 
given by ([2]). This allowed to apply the derivation property of the generator of 
the group 5jFj|j^q to further evaluate its effect on the left factor of A^. 



5 Application to the Wess— Zumino model 

The model we will consider is the simplest Wess-Zumino model, with only one 
type of superfields (and not two as in our preceding work [8j). The general 
non-renormalization properties of such a model ensures that the vertex is not 
renormalized and in fact, in a massless theory, it cannot get any correction. The 
only contribution to the renormalization group comes from the correction to the 
propagator and the supersymmetry has here the effect that all components of 
the superfield get the same renormalization factoi0. We therefore have the 
simple relation 

/3 = 37i. (14) 

so that eq. p3|) becomes 

(A; + l)7fc+i =7i(l + 3aaa)7fe (15) 

^In a supersymmetric gauge theory, the situation would be subtler, since the gauge sym- 
metry and the supersymmetry are not easy to combine: in the component formalism, su- 
persymmetry is not fully explicit and in a superfield formalism, the unconstrained superfield 
describing the gauge supermultiplet is of dimension zero, so any power of this field can appear 
in the renormalized Lagrangian. This can be overcome by an appropriate choice of the gauge 
fixing condition and therefore it does not affect physical properties [13) . 
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The first step in our calculation is to compute the effect of the one-loop 
primitive divergence. There are no primitive divergences at two loops and the 
three-loop one is non-planar and therefore subleading in a large A'^ study. For 
this we will insert the full renormalized propagator at order n in the one loop 
diagram to derive 71 at the order n + 1. Since 7^ is at least of order k, we 
only have to consider a finite expansion in t of the correction factor of the 
propagator. 

We therefore need to compute a simple one loop integral, but with propa- 
gators having logarithmic corrections. All of them can be obtained multiplying 
the propagators r^^(p^) by a factor {p^ / ^i^Y ■ Indeed, differentiating with re- 
spect to X at X = will give the logarithm factors in each of the propagators. 
This procedure is essentially the same as the Mellin transform method used 
in[n]. 

As in [8], we only need to compute the graph for the correction to the aux- 
iliary field propagator. Indeed, the fact that we have corrections on the two 
propagators does not change the observation already made there: the correc- 
tions to the free propagators of all the members of the supermultiplet are equal. 
We are interested in the quantity 



9' f.. 1 



(p2)l-^[(g-p)2]l-y 



= + 8vr^r(l - .)r(i -y) j (^7T^ i"^^ 
= I 9' . 0..+, T{-x-y)T{l + x)T{l + y) 

^87r2 ' T{2 + x + y)T{l - x)T{l - y) ^ ' 

The overall plus sign comes from the product of the minus sign in every cor- 
rection to the IPI two-point function and the free propagator which is —1. It 
will be convenient to introduce the "fine structure constant" a = g'^/Svr^. Now, 
we can get an expression similar (but simpler) to the ones in [11] using the 
expansion of the logarithm of the F function 

ne, y) = -^f^ (2 1 ^(^'"' - - ^''')) <m 

where z = x + y and ({21 + 1) is the Riemann zeta function. From ()17p we 
can obtain the two-point functions calculated with modified propagators by 
applying differential operator in x and y. Differentiation with respect to log(g^) 
produces a factor z which exactly compensates the term which is singular for 
z = 0. Let us define 



H{x,y) 



dT{q^,x, y) 



d\og{q^ 



(18) 



Then, one has 



F(x, y) = 3^^^ exp [2 f; ^|^((x + yr^' - - y''^')\ (19) 
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and 

The evaluation of the 7 function does not require any further regularization. 
After the order a term, which is obtained making x = y = 0, the following terms 
involve corrections to the propagators, which are obtained by taking derivatives 
with respect to x and y evaluated at the origin, this corresponding to diagrams 
with the propagators multiplied by given powers of log(p^). 

We could suppose that all the formalism we developed is superfluous and 
that we could simply derive further eq. (fTTI) with respect to log(g^) to obtain 
directly the higher terms in the expansion. Although this is possible, the com- 
putation are really much more cumbersome, even for the leading term of the 
derivatives of Furthermore, one would obtain the development of the 2- 
point function in the effective action, so that an inversion of the expansion 
would be needed to obtain the propagator. 



6 Calculation and results 

Once the Mellin transform ()19p is known, it is straightforward to obtain the first 
orders in the development of the 7 function from eqs. (jlSp and (j20p . We give 
in tabled] the 12 first terms of the expansion of 7. The principal limitation for 
calculating this "exact" form for each term stems from the rapid growth, which 
behaves as exp(Q;n^/^), of the number of contribution one has to take into ac- 
count. Already in the terms presented in table [H one can see the rapid growth 
of their size. At order 25, there are already 122 possible products of C values 
contributing to the result and this number reaches 409174 at order 100. Such 
a number, as well as the size of the corresponding rational coefficients, limit to 
less than 30 the number of terms which can be computed with a straightfor- 
ward implementation on today's average computer, mainly due to the memory 
footprint. The rapid growth of the required resources limit to a few units the 
number of additional terms which can be computed using more efficient coding 
schemes and bigger memories. The highest orders of the development cannot 
therefore be reached in this way, all the more so that at the end we will need 
to convert these results to some concrete numerical approximation. 

It is therefore convenient to work from the start with numerical values. 
A much higher number of terms is now easily within reach, since for a fixed 
precision, the computation work grows only quartically with the desired number 
of terms. In Table [21 we give the development of 71 up to 16 loops and the 

7 last terms we calculated. The full table, with all results up to 201 loops, is 
available upon request from the authors. 

In the calculation of the expansion of 71, we first have to determine the 
Taylor expansion of the function H{x,y). This is the part of the calculation 
which is numerically most demanding, since these coefficients are the sum of 
a great number of terms which nearly cancel, so that care must be taken that 
the true precision of the result is lesser than the number of digits used in 
the computation. However, something rather remarkable happens: the Taylor 
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7(a) = a - + Ua^ + (^-160 + 16 C(3) j 
+ (2444 - 328 C(3)) + (-45792 + 7056 C(3) + 2016 C(5)) 
+ (1 005 480 - 169 152 C(3) - 70 896 C(5) + 8960 (3)) 
+ (-25 169 760 + 4 509 408 C(3) + 2 199 840 C(5) 
+564480 C(7) - 390400 C^(3))a^ 
+ (705 321 200 - 132 548 640 C(3) - 69 922 848 C(5) 
-29005 632 C(7) + 14 193 504 C^(3) +6397056 C(3)C(5))a^ 
+ (-21 841 420 384 + 4 261 047 424 C(3) + 2 354 993 856 C(5) 
+1 194 909 696 C(7) + ^^^^f C(9) - 512 441 536 C^(3) 
-383 788 416 C(3) C(5) + ^^^^ (^{3)^ a^^ 
+ (740194188032 - 148 784410432 C(3) - 84 779 661 888 C(5) 
-47 818 582 272 C(7) - 58 999 853 440 ^^q^ ^ 19225 297088 C^(3) 
+17 828 697 216 C(3) C(5) + 1 829 076 480 (5) 
+3 838 602 240 C(3)C(7) - 3432 237056 ^8(3)^^11 

+ (-27 243 674 154 368 + 5 610 375 120 768 C(3) + 3 266 192 145 024 C(5) 
+1961976190464 C(7) + 1 019076 124 160 C (9) + 230 546 534400 C(ll) 
-760 702109184 C^(3) - 788057929 728 C(3)C(5) - 141297435 648 C^(5) 
-297887016 960 C(3) C(7) + 59 550068 736 C^(3) 
+34 512 334 848 C^(3) C(5)) a^^ + ■ ■ ■ 

Table 1: The first 12 terms of the development of 71 
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j(a) = a - 2 + 14 - 140.767089549446491434 a"^ 
+2049.72533576365307439 - 35219.8401369368689507 
+741582.310142069315875 - 1.74630317191742523615 x lO'^ 
+4.72719801334671229530 xlO^ - 1.39759545666280992694 x 10^° a^' 
+4.60146704077682933925 xlO^^ a^^ - 1.63220296094286720854 x 10^^ a 
+6.32651854893093835423x101"^ 013- 2.61715263667021333524x10^6 a 
+1.16791189443603376676x10^8 a^^ - 5.52247245848724267096 x lO^^ a 

+ ••• 

+8.4053176185682527526 X 10"^^*^ a^^^ - 4.9339110330514367678 x lO^^'^ 
+2.9112362346747106444x10^6° a^^^ - 1.7263592738217495952x10^^3 
+1.0289894774008300571x10^^66 ai99 _ 6.1636327768018535021 x 10^68 

+3.7107878544109289292 x 10^^^ a^oi + • • • 



Table 2: Approximations of the development of 71 



expansion of the rather complicated function H in (|19p is, asymptotically, the 
same as the one of h, defined by: 

h{x, y) = {l + xy) {-^ + ^ - 1) + 1 (21) 
Vl + x l + y J ^1 - X - y 

More precisely, the terms of total degree n in the Taylor expansion of h and the 
original function H are the same up to a relative error which is smaller than 
2~". If it is relatively easy to understand the properties of the terms in xy"' 
from the serie with general term + 1) — 1, the interplay of different terms 
which yields such a simple result for the general terms appears quite miraculous. 
Even more, the error term proportional to 2~", 3~" and 4~" can be given similar 
explicit forms and the determination of error terms proportional to powers of 
bigger integer seem only to require more terms with higher precision. These 
asymptotic formulas have been used to obtain the higher degrees of the Taylor 
expansion for our numerical application. 

The asymptotic behavior of the coefficients gn in the expansion of 71 in 
powers of a can be determined. With 

00 

71(a) = ^ff„a", (22) 

n=l 

we have that 

g„a (-3)"r(n + 2/3). (23) 

This implies that a suitably defined Borel transform will have a single pole at 
—1/3 as dominant singularity. However the asymmetry between odd and even 
terms of the series implies that there is also a weaker singularity at +1/3. It is 
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not a simple pole but the singularity of the dilogarithm Li2. This precludes the 
analytic continuation of the Borel transform on the positive real axis and its 
approximation by Pade methods. A brutal answer to this problem is simply to 
integrate the Borel transform only on the interval [0, |]. This gives an analytic 
function which has the desired asymptotic expansion at the origin. However 
this procedure is rather ad hoc and gives a far from unique solution. 

There is also the possibility that higher order corrections to the Schwinger- 
Dyson equation stemming from the addition of multiloop primitive divergences 
would weaken the singularity on the positive axis, so that the singularity on the 
positive axis would disappear in a complete calculation. 

Taken as an asymptotic formula, eq. (|2ip could be useful for a demonstration 
of the asympotic properties of the expansion of 71. The second term in h has 
positive Taylor coefficients. Combined with the alternating signs of the main 
term in the development of 71, this introduces strong compensation between 
the terms with even total powers of derivatives in eq. ()20p and those with odd 
total powers. As a result, in the calculation of the term of degree n, the terms 
with one of the propagator with the correction of degree n — 1 and the other 
one equal to the free propagator are dominant. 

7 Discussion 

High orders of the perturbative development of the renormalization group func- 
tion have been obtained, resumming all diagrams with the simplest primitive 
divergence. The number of terms calculated in the corresponding expansion 
allows to reach a clearly asymptotic regime. The Borel transform of the result- 
ing series however presents singularities on the positive real axis. The question 
remains open whether this signals a fundamental limitation of the perturbative 
calculations or it is just an artifact of the present approximation. Indeed, in a 
large N approximation, other primitive divergences are possible, beginning at 
the fourth order. 

However, one should first test if the additional contributions remain sub- 
dominant. Indeed, in the contribution from the four loop primitive divergence, 
the corrections to the propagators can be split between the eleven propaga- 
tors, generating a rapidly growing number of terms and this could more than 
compensate for the fact that corrections to the propagator of order n — 3 are 
of order with respect to the propagator of order n which gives the main 

contribution in the term we have calculated. This would also be an important 
step in the way of conjecturing ("-function style functional equations generaliz- 
ing the one obtained in [11] for the simple case where the 7 function satisfies a 
differential equation. 

One interesting feature of the approach we presented is that it provides 
a clear path for evaluating the correction stemming from additional terms in 
the Schwinger-Dyson equation. The evaluation of the primitive diagrams looks 
formally as the analytic regularization proposed long ago [HI [151 US] but the 
fact that we need such an evaluation only for primitively divergent diagrams 
avoids much of the technicalities stemming from the complex multidimensional 



10 



pole structures appearing in general diagrams. 
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